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Intersections of thick, plane SU{2) center vortices are characterized by the topological 
charge \Q\ — 1/2. We compare such intersections with the distribution of zeromodes of the 
Dirac operator in the fundamental and adjoint representation using both the overlap and 
asqtad staggered fermion formulations in SU (2) lattice gauge theory. We analyze configura- 
tions with four intersections and find that the probability density distribution of fundamental 
zeromodes in the intersection plane differs significantly from the one obtained analytically 
in [ij. The Dirac eigenmodes are clearly sensitive to the traces of the Polyakov (Wilson) 
lines and do not exactly locate topological charge contributions. Although, the adjoint Dirac 
operator is able to produce zeromodes for configurations with topological charge \Q\ = 1/2, 
they do not locate single vortex intersections, as we prove by forming arbitrary linear com- 
binations of these zeromodes - their scalar density peaks at least at two intersection points. 
With pairs of thin and thick vortices we realize a situation similar to configurations with 
topological charge \Q\ = 1/2. For such configurations the zeromodes do not localize in 
the regions of fractional topological charge contribution but spread over the whole lattice, 
avoiding regions of negative traces of Polyakov lines. This sensitivity to Polyakov lines we 
also confirm for single vortex-pairs, i.e., configurations with nontrivial Polyakov loops but 
without topological charge. 
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I. INTRODUCTION 

Quantum chromodynamics (QCD) at low energies is dominated by the non-perturbative phe- 
nomena of quark confinement and spontaneous chiral symmetry breaking (SCSB). Presently, a 
rigorous treatment of them is only possible in a lattice regularization. Many of the important 
features of non-abelian gauge theories are already present in SU{2), which simplifies theoretical 
and numerical calculations. 

The non-perturbative vacuum can be characterized by various kinds of topological gauge field 



excitations. An established theory of SCSB relies on instantons 



m 



51] , which are localized in space^ 



time and carry a topological charge of modulus one. Due to the Atiyah-Singer index theorem 



a 



0- 



a zeromode of the Dirac operator arises, which is concentrated at the instanton core. In the 
instanton liquid model 



101, 



111 ] overlapping would-be zeromodes split into low lying non-zeromodes 



which create the chiral condensate. 

On the other hand there is plenty of evidence 



12 



17l | for the explanation of confinement by 



center vortices, closed, quantized magnetic fiux tubes with values in the center of the gauge group. 



These properties are the key ingredients in the vortex model of confinement, w. 
appealing and was also confirmed by a multitude of numerical calculations [l8 



lich is theoretically 
Lattice simu- 



20|-l23|], 



lations indicate that vortices may be responsible for topological charge and SCSB as well 
and thus unify all non-perturbative phenomena in a common framework. However, to date the 
potential physical mechanism for symmetry breaking through vortices is still unclear. A similar 
picture to the instanton liquid model exists insofar as lumps of topological charge arise at the 
intersection and writhing points of vortices [24. |25|. 



The present numerical investigation concentrates on the topological charge contributions of vor- 
tex intersections and their localization by Dirac zeromodes. All measurements were performed on 
hyper-cubic lattices of even sizes from 12"^ up to 22'^-lattices. We start with the construction of thick 
planar vortex configurations and derive the topological charge contribution of their intersections. 
Then we introduce the lattice index theorem for various fermion realizations and representations 
in more detail. 

Using the overlap and asqtad staggered Dirac operator, we compute fundamental zeromodes in 
the background of four vortex intersections. By visualizing the probability density, we compare 
the distribution of the eigenmode density with the position of the vortices and the topological 
charge density created by intersection points. For interpreting the results, we refer to existing 
analytical calculations of the zeromodes for fiat vortices with gauge potentials living in the Cartan 



subalgebra of a SU{2) gauge group A comparison reveals systematic discrepancies, which we 
tentatively attribute to differences in the values of the Polyakov loops between the analytical [l| 
and our numerical calculations. We argue that configurations with the same field strength, but 
different Polyakov loops do not give rise to the same eigenmode density. Comments on the physical 
significance of this observation follow further below. 

Furthermore, we calculate the Dirac eigenmodes in the adjoint fermion representation. Adjoint 
fermions are of special interest with respect to the constituents of the QCD vacuum with fractional 
topological charge. For configurations with total topological charge \Q\ = 1/2 no fundamental 
zeromode would necessarily be produced, however, adjoint fermions are able to create a zeromode. 



Edwards et al. presented in 
Garcia-Perez et al 



26( 1 some evidence for fractional topological charge on the lattice. 
27l |. however, associated this to lattice artifacts, i.e., topological objects of size 
of the order of the lattice spacing. 

Vortex intersections are examples of fractional topological charge contributions \Q\ = 1/2, which 



can be related to merons [281] and calorons [29(]. Nevertheless, on lattice configurations with periodic 
(untwisted) boundary conditions, no single vortex intersection can be realized. Therefore we tried 
to separate adjoint zeromodes on periodic lattices with four vortex intersections in order to find 
linear combinations of zeromodes which detect one vortex intersection only. We did not find such 
linear combinations, the scalar density of the zeromodes peaks at least at two intersection points. 

Finally we combine "thin" and "thick" vortex sheets, which simulate in some sense twisted 
boundary conditions. In fact, "thin" vortices are not recognized by adjoint fermions. If we intersect 
such vortex pairs adjoint fermions detect a single vortex intersection only and the adjoint index 
theorem signals the topological charge Q = 1/2. Nevertheless, the zeromodes do not seem to 
localize the vortex intersection but rather extend over the whole lattice, avoiding regions with 
nonvanishing topological charge density and negative Polyakov lines. 



II. PLANE VORTICES 



We investigate planar vortices parallel to two of the coordinate axes in SU{2) lattice gauge 
theory. Since we use periodic (untwisted) boundary conditions for the links, vortices occur in pairs 
of parallel sheets, each of which is closed by virtue of the lattice periodicity. We use two different 
orientations of vortex sheets, xy- and zt-planes with nontrivial links within the vortex thickness. 
These links vary in the a-s subgroup of SU{2), Ufj_ = exp{i(j)a3). For xy-vortices fj, = t links are 
nontrivial in one i-slice only, for zt-vortices we have nontrivial y-links in one y-slice. Since the 
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U{1) subgroup remains unchanged, the direction of the flux and the orientation of the vortex are 
determined by the gradient of the angle (p, which we choose as a Hnear function of the coordinate 
perpendicular to the vortex. Upon traversing a vortex sheet, the angle <j) increases or decreases by 
vr within a finite thickness 2d of the vortex, see Fig.[TJ Center projection leads to a (thin) P-vortex 
at half the thickness (d) [18|. We distinguish parallel and anti-parallel vortices, i.e., vortex sheets 
with the same resp. opposite orientation. The angle <j)i first increases from to vr, then (pi increases 
to 27r and <j)2 returns to 0. For an xy- vortex with vortex sheets at zi and Z2 the t-links in one t-slice 
vary with the angle 
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1 {4'2) points in the same (opposite) direction at the two vortex sheets of a 



Since the gradient of 

pair, their fiuxes are (anti-) parallel and the total flux through the zt-plane is therefore 2vr (zero) 




FIG. 1. The link angle a) 4>i of a parallel and b) 02 of an anti-parallel vortex pair. The arrows rotate 
counterclockwise with increasing 0^. The vertical dashed lines indicate the positions of the P- vortices. In 
the shaded areas the links have positive, otherwise negative trace. 



Next we consider these thick, planar vortices intersecting orthogonally. As shown in [30|], each 
intersection carries a topological charge with modulus \Q\ = 1/2, whose sign depends on the 
relative orientation of the vortex fluxes. The plaquette definition simply discretizes the continuum 
(Minkowski) expression of the Pontryagin index to a lattice (Euclidean) version of the topological 
charge definition: 



Q 



1 



16vr2 



■-^^ I d'^xef,uai3tT[J='ai3J^iJ.u] = ^ j d'^xE-B (2.2) 
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We can drop color indices since all our links belong to the as-generated U{1) subgroup of SU{2). 
Our xy-vortices have only nontrivial zt-plaquettes, i.e., an electric field Ez, while zi- vortices bear 
nontrivial xy-plaquettes corresponding to a magnetic field B^- The topological charge is then 
proportional to E^B^, hence parallel crossings give Q = 1/2 and anti-parallel crossings give Q = 
— 1/2. However, due to the finite lattice spacing (or the finite thickness of vortices in units of the 
lattice spacing) the numerical absolute values of the lattice charge are slightly smaller than the 
continuum ones. 

As mentioned above, the lattice periodicity forbids single vortex sheets and therefore we get 
at least four intersections for two vortex pairs, summing up to an even valued topological charge. 
We combine xy and zt-vortices in central y- and i-slices with vortex centers at xi^2 resp. zi^2 
located symmetrically around the lattice center and varying vortex thickness d. In Fig. [2] we 
present a 3-dimensional view of the intersecting vortices on a 12^-lattice together with topological 
charge distributions in the xz-plane at (y = 6,t = 6), which is the intersection plane. Further 
density plots of Dirac eigenmodes or topological charge, as well as Polyakov loop distributions 
will all be plotted in this plane, except for (specially mentioned) orthogonal plots in order to 
analyze localization properties. The four intersection points of two parallel vortices all carry 
topological charge Q = +1/2 whereas for two anti-parallel vortices and parallel-anti-parallel vortex 
combinations we get two intersections with Q = +1/2 and two with Q = —1/2. 



Parallel Vortices 



Geometry 



Anti-parallel Vortices 




FIG. 2. A S-dimensional section (hyperplane) in xyz-direction of a 12^-lattice at time t = 6 (center). The 
horizontal planes are the xy-vortices, which exist only at this time. The vertical lines are the zt-vortices, 
which continue over the whole time axis. The ticks protruding from the vertical lines extend in time direction. 
The vortices intersect in four points of the y = t = Q - plane, giving topological charge Q = 2 for parallel 
vortices (Ihs) or Q = for anti-parallel vortices (rhs). 
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III. FERMIONIC ZEROMODES OF THE OVERLAP AND ASQTAD STAGGERED 
DIRAC OPERATOR FOR INTERSECTING CENTER VORTEX FIELDS 

We test the lattice index theorem and analyze the scalar density p{x) = i/j^'djh:) of fermionic 



31[ the improved 



zeromodes ip in the background of intersecting plane vortices. As described in 
staggered operator also produces eigenmodes which can clearly be identified as zeromodes and all 
results in this paper show perfect agreement between the two fermion realizations. The fermionic 
zeromodes are used to measure the topological charge Q via the Atiyah-Singer index theorem 

m 

ind D[A]=n_-n+ = Q, (3.1) 

where n_ and are the number of left- and right-handed zeromodes of the Dirac operator D. 
This equation accounts for Wilson and overlap fermions in the fundamental representation. The 
adjoint version of the index theorem reads 

ind D[A] = - n+ = 2NQ = 4Q (3.2) 

where = 2 is the number of colors and the additional factor 2 is due to the fact that the fermion is 
in the real representation, hence the spectrum of the adjoint Dirac operator iD is doubly degenerate. 
The eigenvalues of the staggered fermion operator have a twofold degeneracy due to a global charge 
conjugation symmetry in SU{2). We therefore have ind D[A] = n__ — n+ = 2Q for fundamental 
and ind D[A] = n_ — ?i+ = 8Q for adjoint (asqtad) staggered fermions. 



A. Fundamental zeromodes for intersecting center vortex fields with topological charge 

g = 2 

We intersect two parallel vortex pairs with xi = zi = 6 and X2 = ^2 = 16 at y = t = 11 
respectively on a 22^-lattice. The four intersection points all carry topological charge contributions 
of +1/2 and therefore sum up to a total topological charge Q = 2. In agreement with the lattice 
index theorem we get two overlap and four asqtad staggered zeromodes of negative chirality (left- 
handed) in the fundamental representation. Fig.[3]shows the scalar density plots of the fundamental 
overlap and asqtad staggered zeromodes with periodic boundary conditions together with the sum 
of Wilson lines in y- and t-direction (Polyakov-loops) in the intersection plane as well as the scalar 
density plot of the two overlap zeromodes with usual antiperiodic boundary conditions in time 
direction (asqtad staggered modes again distribute similarly). The individual modes all distribute 
equally, showing four distinct maxima, as trivially all their linear combinations do. A close look 
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shows that the zeromodes do not exactly peak at the vortex intersections, they rather avoid regions 
with negative Polyakov hnes and approach the intersections (or the vortex surfaces) from regions 



32] for spherical vortices, 



with positive Polyakov lines. This behavior was already observed in 
the zeromodes avoid regions with negative Polyakov lines. The sensitivity of Dirac eigenmodes 
to the Polyakov value can be exhibited clearly using one parallel vortex pair (e.g. two parallel 
xy- vortices), which are exchanged by a discrete translation in the z-direction by half the lattice 
length Nz- The configuration apparently consists of two identical flux lines because the field 
strength is invariant under the z-translation T^. However, this operation is not a symmetry of 
the Dirac operator since it would change the Polyakov loop, which is a gauge-invariant quantity. 
Therefore the distribution of the scalar density of the Dirac eigenfunctions is not invariant under 
Tz- In fact, the antiperiodic boundary conditions in time direction change the time-like Polyakov 
lines and the zeromode density shows the shift in the z-direction by half the lattice size, see Fig. [3ji. 

The configuration is also presumably equivalent to the abelian gauge fields whose zeromodes 
are derived analytically in Nevertheless, the probability distribution in the intersection plane 
shown in the reference differs significantly from the one we obtained, namely in that the former 
exhibits exactly the translation symmetry discussed above. This may indicate a relation to the issue 
of the Polyakov loop because configurations having the same flux distribution may still differ in the 
values of their Polyakov loops. The values of the Polyakov loops in our study and the calculation of 
[l| only coincide in the intersection plane. The choice of boundary conditions is another possibility 
for the origin of the discrepancy. For the flux of two parallel center vortices the transition matrices 
need to have a winding number one around the boundary. On the lattice, the transition functions 
are supposed to be incorporated in the links such that all plaquettes agree with the corresponding 
integrals of the continuum field. Hence it should be sufficient to use periodic boundary conditions 
for the links as well as for the fermions. However, this allows further zeromodes which cancel each 
other in the index theorem (non-topological zeromodes) and so we prefer the usual antiperiodic 
boundaries in time direction which only reproduce the relevant zeromodes. 



B. Adjoint zeromodes for intersecting center vortex fields with topological charge Q = 

and g = 2 

Now we try to locate the fractional topological charge contributions with adjoint eigenmodes of 
the overlap Dirac operator which are also sensitive to topological charge contributions of \Q\ = 1/2. 
For asqtad staggered fermions we find the correct (doubled) numbers of zeromodes in all cases and 
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FIG. 3. Q = 2 configuration: Scalar density plots of a) the two overlap zeromodes, b) the four asqtad 
staggered zeromodes, both with periodic boundary conditions, and d) the two overlap zeromodes with 
antiperiodic boundary conditions in time direction. The plot titles indicate the plane positions, the chirality 
(chi=±l) and numbers (n=l-2/n=l-4) of (overlap/asqtad staggered) zeromodes and the maximum density 
peak in the plot, c) Sum of Wilson lines in y- and t-direction (Polyakov-loops) in the intersection plane. 
P-vortices are indicated with red or black lines. 



the scalar densities of the sum of all zeromodes always localizes similarly. 

We do not find zeromodes localized to a single region with nonvanis. 
contributions. Therefore we use the inverse participation ratio (IPR) 



ring topological charge 



3314361] to quantify the 



localization of eigenmodes. The IPR of a normalized i^^Pii^) ~ ^) fi^ld Pi{x) is defined as 

TV 



I = N^pUx) 



(3.3) 



x=0 



where N is the number of lattice sites x. With this definition, / characterizes the inverse fraction of 
sites contributing significantly to the support of p{x), i.e., a high IPR indicates that the eigenmode 
is localized to a a few lattice points only. We perform systematic and random IPR maximization 
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procedures for linear combinations of zeromodes in order to get single eigenmode peaks localized 
to regions with nonvanishing topological charge contribution. 
Q=0: 

For the configuration with topological charge Q = we intersect two anti-parallel vortices 
with the same vortex centers as for the Q = 2 configuration described above (xi = zi = 6 resp. 
X2 = Z2 = 16) at y = t = 11 respectively on a 22^-lattice. For this configuration we do not get 
fundamental zeromodes, but with periodic boundary conditions we find two adjoint overlap (four 
asqtad staggered) zeromodes of each chirality. Fig. H] shows the scalar density plots of the adjoint 
overlap zeromodes for the Q = configuration. The left-handed zeromodes (Fig. H^) peak at the 
intersection point (6, 11, 6, 11) of topological charge Q = 1/2 with a maximum density of 8.595-10"^ 
but a second maximum can be found near the intersection plane at (16, 12, 6, 12) with a maximum 
value of 8.601 • 10~® (Fig. HJa) which is next to the intersection with Q = —1/2 at (16, 11,6, 11). 
These zeromodes rather localize the xy- vortex sheet at 2 = 6, pronouncing the intersections whereas 
the right-handed zeromodes of Fig. locate the other xy- vortex sheet at z = 16, but with notches 
at the two intersection points. The individual zeromodes of the same chirality show identical scalar 
density distributions, therefore it is not possible to find linear combinations of the zeromodes which 
locate single vortex intersections. 




FIG. 4. Q = configuration: Scalar density plots of the two adjoint overlap zeromodes (n=l-2) with negative 
chirality (chi=-l) in the .xy-plane at a) y = t = 11 (intersection plane, red lines indicate the P-vortices), b) 
y = t = 12 and c) the two right-handed (chi=l) adjoint overlap zeromodes at y = t = 11. 

Q=2: 

The configuration with topological charge Q = 2 on the 22^-lattice from above gives 16 adjoint 
asqtad staggered zeromodes with negative chirality for antiperiodic boundary conditions in time 
direction. For adjoint overlap fermions we use a 16^-lattice with vortex centers xi = zi = A and 
3^2 = -22 = 12 at y = t = 8 respectively and we find eight left-handed zeromodes. With periodic 
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boundary conditions we find two positive (non-topological) and ten negative zeromodes. Fig. [5] 
shows the scalar density plots of adjoint overlap and asqtad staggered zeromodes of the Q = 2 
configuration. The zeromodes locate the xy-vortex pairs but even individual modes do not show 
single peaks locating one intersection with Q = 1/2. Linear combinations of the eight zeromodes 
with negative chirality show six distinct IPR maxima. This was obtained by a systematic study 
with the eight coefficients of the linear combination varying from —10 to 10 in integer steps. 
Further we started from 20.000 random points in the parameter space and determined the nearest 
maximum by the gradient method. Each of the six maxima was found between 2.000 — 6.000 
times and no other maxima were obtained. In Fig. [5^ we plot a 2D cut through the first three 
IPR maxima in the 8D parameter space of linear combinations. The scalar density of the linear 
combination of the eight zeromodes with maximal IPR is presented in Fig. [Sf, it still peaks at two 
vortex intersections. The ten zeromodes for periodic boundary conditions look more symmetric, 
but again no single mode locates one vortex intersection. We again performed a systematic search 
with the ten coefficients of the zeromodes varying from —6 to 6 in integer steps and again 20.000 
random points. In this case we found only five maxima, each 3.000 — 5.000 times. Again, it is not 
possible to find linear combinations of the zeromodes to locate at single vortex intersections. 



C. Adjoint zeromodes for center vortex fields with one single "thick" intersection and 

topological charge (5 = 1/2 

We analyze a configuration of "thin-thick" vortex intersections apparently having topological 
charge \Q\ = 1/2. The profile of a "thin-thick" vortex is plotted in Fig. [6^. The thin vortex sheet 
is defined by the jump of the y- or t-link from -|-1 to —1 at the boundary. The thick vortex is 
located symmetrically around the center of the lattice with thickness d. The thin-thick xy- and 
zt-vortices still intersect at four points, but the plaquette or hypercube definitions of topological 
charge do not recognize the thin vortex sheets and therefore only measure one topological charge 
contribution Q = 1/2 of the "thick" vortex intersection, see Fig. [6}). 

We compute the Dirac eigenmodes on a 22^-lattice with two intersecting " thin-thick" xy- and zt- 
vortices with vortex thickness d = 20att = y = ll, respectively. We get no fundamental zeromodes 
for the above topological charge "|(3| = 1/2 configuration". But the adjoint Dirac operator does 
not recognize the thin vortices, just like the field theoretic operators we used to measure topological 
charge density. We find two adjoint overlap and four adjoint staggered zeromodes with negative 
chirality, which due to the index theorem again results in Q = 1/2. For the adjoint fermions this 
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y=ll, t=ll, chi=-l, n=l-16, max=0.001 15675 y=8, t=8, chi=-l, n=l-8, max=0.00 196495 








16 16 



e) f) 

FIG. 5. (5 = 2 configuration: Scalar density plots of the a) 16 asqtad staggered and b) eight overlap 
left-handed adjoint zcroniodes for antiperiodic boundary conditions, the modes locate the xy-vortex pair. 
Further, c) two right-handed (non-topological) and d) ten left-handed adjoint overlap zeromodes for periodic 
boundary conditions, e) Plane through the highest three IPR maxima in the parameter space of linear 
combinations of the eight zeromodes. The peaks are very broad and easy to identify, f) The scalar density 
of the linear combination of the eight zeromodes with maximal IPR still peaks at two vortex intersections. 
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FIG. 6. "Q = 1/2 configuration": a) Link profile of a "thin-thick" plane vortex, the link angle (f) (blue) 
decreases from tt to within a certain vortex thickness d. The arrows (links) rotate counterclockwise with 
decreasing 0. The thin vortex is given by the jump at the boundary. The red dashed line shows the trace 
of adjoint links TvUa (see text below) b) Topological charge density in the intersection plane. 

configuration truncates three of tlie four vortex intersections and in this way simulates a situation 
related to the one achieved by twisted boundary conditions, namely, a single detectable intersection. 
Therefore it is possible to have a configuration that looks like having fractional topological charge. 
The eigenmode density distributions of overlap and asqtad staggered zeromodes are identical, Fig. [7] 
shows the former for the intersection plane and for orthogonal planes to it at the intersection point. 
The modes are clearly sensitive to the traces of the adjoint link Ua (see also Fig. [6^) 

(TrC/)(TrC/)t = {TrUf = 1 + TtUa, (3.4) 

which in our case define the adjoint Polyakov lines Pa (Wilson lines), since we have only nontrivial 
t-links in one time-slice and nontrivial y-links in one y-slice. The zeromodes prefer regions of 
positive Polyakov lines and avoid negative Polyakov lines with respect to the boundary conditions. 
For the antiperiodic boundary conditions in time direction the signs of the Polyakov lines are 
exchanged. Hence, the zeromode densities peak at the xy- vortex center at z = t = 11 and avoid 
the zt- vortex center at x = y = 11, or rather peak at the boundary parallel to the zt- vortex at 
x = or a; = 22 in the y = 11-plane. For completeness we should mention that we get exactly the 
same results if we analyze a "Q = —1/2 configuration", which can be realized if we rotate the links 
of one of the thick vortices in the opposite direction. In this case we get two adjoint overlap and 
four adjoint staggered zeromodes of positive chirality which show the same behavior. 

We further vary the vortex thickness d, see Fig. [8] and do find analogue results: The adjoint 
zeromodes approach the thick vortex intersection from regions of positive Polyakov lines, but do 
not localize exactly at the region with topological charge contribution \Q\ = 1/2. They rather 
spread over the whole lattice avoiding regions of negative traces of adjoint Wilson lines. 
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y=ll, t=ll, chi=-l, n=l-2, max=0.000070289 "=11. ^=11' chi=-l, n=l-2, max=0.0000192506 





FIG. 7. "Q = 1/2 configuration": Scalar eigenmode density of two adjoint overlap (identical to four 
adjoint asqtad staggered) zeromodes of negative chirality in various planes through the thick-thick vortex 
intersection: a) xz-plane (intersection plane): The zeromodes avoid regions of negative adjoint Polyakov 
(Wilson) lines {Pa, red dots) with respect to boundary conditions (see text above, Eq. (|3.4p and Fig. |6^) 
and therefore do not peak at the topological charge contribution Q — 1/2 of the "thick-thick" vortex 
intersection at a; = z = 11; the "thin-thin" and "thin-thick" intersections are not recognized by the Dirac 
modes; b) j/i-plane: The antiperiodic boundary conditions invert the profile in time- compared to the one 
in y-direction, hence the zeromodes prefer the xy- vortex at t = 11 and avoid the zt- vortex at y = 11 (the 
lines indicate the xy (red) and zt (green) vortex sheets in the i = 11 resp. y = 11 slices); c) xy-plane: The 
zeromodes reflect the profile of the adjoint y- Wilson lines {Pa, red dots) of the zi- vortex in cc-direction at 
y ~ 11; d) zt-plane: The zeromodes refiect the profile of the adjoint Polyakov lines {Pa, red dots) of the 
a;2/- vortex in z-direction at t = 11, inverted by the antiperiodic boundary conditions in time direction. 
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y=ll, t=ll, chi=-l, n=l-2, max=0.0000706566 y=ll, t=ll, chi=-l, n=l-2, max=0.000084769 




y=ll, t=ll, chi=-l, n=l-2, max=0.000107012 y=ll. t=ll. chi=-l, n=l-2, max=0.OO0163089 




FIG. 8. "Q = 1/2 configuration": Scalar cigcnmodc density of two adjoint overlap (n=l-2) or four adjoint 
asqtad staggered (n=l-4) zeromodes of negative chirality in the intersection plane of zt- and xy- vortices for 
various vortex thicknesses: a) dx ^ dz = 16 b) dx = dz = 12 c) dx = dz ^ 8 d) dx = dz = A d) dx = 20, 
dz ~ 8 f) dx — 8, dz = 2. The adjoint zeromodes approach the thick vortex intersection at x = 2 = 11 
from regions of positive, adjoint Polyakov lines (Pa), but do not localize exactly the topological charge 
contribution (5 = 1/2 because they strictly avoid regions of negative, adjoint Polyakov (Wilson) lines. We 
also plot the profile of the adjoint Polyakov (Wilson) lines {Pa, red dots). 
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IV. DIRAC MODES IN THE BACKGROUND OF SINGLE CENTER VORTEX PAIRS 

Since we found that the Dh'ac zeromodes seem to be more sensitive to the Polyakov (Wilson) 
hnes than to topological charge contributions we analyze configurations of single vortex pairs 
apparently without topological charge. We use parallel xy-vortices with t-links in one time-slice 
(to) varying in z-direction according to Eg. 12.11 shown in Fig.[T^). Due to the translation symmetry 
in X and y, i.e., parallel to the vortex planes, and the fact that the links vary only within an abelian 
U(l)-subgroup of SU{2), our configuration corresponds to an abelian 2-dimensional problem which 
we can compare wi_th the abelian center vortices on considered in {l|. The nontrivial transition 



functions used in [l| are incorporated in our nontrivial links and the analytical result should be 
equivalent to our configuration with periodic fermion boundary conditions. 

The occurrence of zeromodes can be related to a U(l) index theorem for T^, Ref. [l| finds for 
m+ vortices and m_ anti-vortices a number Am/2 = Irri-L. — m_|/2 of zeromodes (for more details 

n 

see also [37|]). Since we work in four dimensions and with SU{2), we expect to find four times 
more zeromodes: we have two times more color indices and two times more spinor components. 
This is indeed our result for periodic boundary conditions in the x- or y-directions and antiperiodic 
boundary conditions in the z- or t-directions, i.e., perpendicular to the vortex surfaces Q. 

We find two overlap zeromodes of each chiralitjo spreading over the whole lattice parallel to the 
vortex surfaces, i.e., their scalar density distribution is constant in the x- and y-direction, strictly 
avoiding regions with negative Polyakov values. At the time-slice with nontrivial t-links (to) the 
density distribution clearly follows the profile of these links, which correspond to the Polyakov 
loops. We confirm these results also for staggered fermions and adjoint representations with ad- 
joint Polyakov traces Pa = P"^ - I (see Eq. [MD- Fig. [9] shows our Polyakov profiles P and Pa 
together with scalar density plots of the fundamental and adjoint Dirac zeromodes with periodic and 
antiperiodic boundary conditions in t-direction (the latter invert the Polyakov profiles, correspond- 
ing to a multiplication with —1). The individual zeromodes for both chiralities show equivalent 
density distributions, but they differ from the analytic results found in [l[. Further, the lowest 
non-zeromodes also propagate parallel to the vortex surfaces, showing similar density distributions 
but with opposite response to the Polyakov profiles, i.e., exchanged boundary conditions. 

The above configuration consists of two identical fiux lines because the field strength is invariant 



^ For antiperiodic boundary conditions in the x- or y-directions, i.e., parallel to the vortex surfaces we do not find 

zeromodes. Because of the translation invariance in x and y, the Dirac modes should contain a free wave factor 

with a wave vector k = {kx,ky). But for antiperiodic boundary conditions the smallest allowed modulus of k is 

^ 7^ and the corresponding lowest eigenvalue must be greater than zero. 
^ With the introduction of a little noise, these zeromodes would mix and be lifted to low lying modes. 
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X=6, y=6, chi=-l, n=l-2, max=0.000221825 




y=6, fee, chi=-l, n=l-4, max=0.000267027 



x=6, y=6, chi=-l, n=l-4, max=0.000683854 




0.0006 
0.0004 
0.0002 




FIG. 9. Scalar density plots of Dirac zcroniodes in the background of a single xy- vortex pair at t = 6 
and zi^2 = 3.5,9.5 (see Eg. 12. ip : a,c) two fundamental zeromodes for each chirality (the four individual 
modes have equivalent density distributions) show the Polyakov profile at t = 6 and arc constant in x- (y-) 
direction; b) profile of the (adjoint) Polyakov traces P (Pa) for our configuration (indicated with red dots in 
the density plots); d) two fundamental zeromodes for each chirality with antiperiodic boundary conditions 
show the inverted Polyakov profile att~6; e,f) four adjoint zeromodes for each chirality (the eight individual 
modes again have equivalent density distributions) with periodic resp. antiperiodic boundary conditions. 
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under a z-translation Tz by half the lattice size, which exchanges the two vortices. But this 
symmetry is spoilt by the fact that in the zt-plane we need to distinguish two regions which differ 
by the value of the Polyakov loop P. More generally, we could multiply all links Ut at a fixed 
t-coordinate (not necessarily to) by an arbitrary 5'C/(2)-element. While this would not alter any 
plaquettes (i.e., the flux), P would change at a given point z and could be modified to any desired 
value. Invariably, however, the Polyakov loop varies by a center element at the locations of the 
vortices. This simply follows from Stokes theorem. Two Polyakov loops on either side of the vortex 
can be deformed into a single curve encircling it. The line integral along this curve equals the center 
element of the vortex in its interior. 

Consequently, corresponds to a center transformation switching the sign of P. Since the 
Polyakov loop is a gauge-invariant quantity, the configurations before and after a translation with 
Tz cannot be gauge-equivalent and in general the Dirac operator will not be invariant under Tz- 
Since the scalar densities of zeromodes and the near-zeromodes are asymmetric with respect to the 
translation Tz and the only effect of Tz is a change of P, we are forced to conclude that the Dirac 
eigenmodes are sensitive to the value of the Polyakov loop. 

For completeness we mention the anti-parallel vortex (Fig. [Tb). Since the total flux vanishes, 
in accord with the 2-dimensional U(l) index theorem [37| we observe no zeromodes. In this case 
the translation Tz changes both, the direction of fluxes and the value of P. Still, the lowest-lying 
eigenmodes are sensitive to the value of the Polyakov loop with respect to boundary conditions. 



V. CONCLUSIONS 



In this paper we have studied low lying asqtad staggered and overlap Dirac operator eigenmodes 
in the background of intersecting vortex configurations. Each intersection region carries topological 
charge ±1/2. We found that the zeromodes are sensitive to the value of the Polyakov loops (Wilson 
lines), avoiding regions with negative Polyakov loops. In particular, we studied configurations 
with two parallel planar vortices, which can be exchanged by a discrete translation T by half the 
lattice size. While the field strength (and the topological charge distribution obtained from the 
field strength) is invariant under this translation, the Dirac operator is not, since the translation 
modifies (shifts) the Polyakov loop. Therefore, the distribution of the scalar density of the Dirac 
operator zeromodes is not forced to be invariant under this translation, and we do find that they 
are not. This is in contrast to the analytical results of Ref. [l| which appear to be invariant under 
the discrete translations. This happens even though our planar vortex configurations are presumed 
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to be equivalent to the abelian gauge configurations studied there. We think the difference might 
be related to differences in the Polyakov loop. Configurations with the same flux distribution can 
differ in the distribution of the Polyakov loop values, which can lead to differences in the scalar 
density of Dirac operator zeromodes, as those tend to avoid regions of negative Polyakov loops. We 
find that the zeromodes in the fundamental representation are not quite spread along the center 
vortex sheets and don't quite peak at the vortex intersection points. Rather, they approach these 
topological structures from regions with positive Polyakov lines while avoiding regions of negative 
Polyakov lines. 

With adjoint fermions we tried to find zeromodes which identify exactly one topological charge 
contribution Q = 1/2 of a single vortex intersection. We analyzed linear combinations of zeromodes 
which maximize the inverse participation ratio (IPR), i.e., localize as much as possible. We found 
that the scalar eigenmode density of adjoint fermions always peaks at least at two intersections. 
Further, we analyzed a lattice configuration with only one "thick" vortex intersection. Since both, 
gluonic and (adjoint) fermionic definitions of topological charge fail to detect intersections with at 
least one thin vortex, both definitions of Q merely signal the value Q = 1/2. We find it remarkable 
that the two corresponding adjoint zeromodes are not localized to the region with nonvanishing 
topological charge contribution but spread over the whole lattice, avoiding regions of negative 
traces of adjoint Polyakov (Wilson) lines. Therefore we conclude that the Dirac zeromodes are 
more sensitive to the Polyakov (Wilson) lines than to the topological charge contributions. 

Finally we confirm this conjecture with Dirac eigenmodes for configurations with nontrivial 
Polyakov profiles of single center vortex pairs, apparently without topological charge. In principle, 
the sensitivity of the Dirac operator to the Polyakov (Wilson) lines with respect to boundary 
conditions is nothing new. _In calorons for example the zeromode hops with the boundary conditions 
in the compact direction 
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39[ between the monopole constituents, where the Polyakov loop 



actually passes through 1 and —1. Since calorons seem only to be important in the deconfinement 
phase, i.e. with finite temporal lattice size, these effects may be interpreted as finite size effects. 
Similarly, our specially constructed vortex configurations discussed in this paper, lacking a natural 
size for a correlation length, could be viewed as having large finite size effects. Thus our observations 
may not completely hold in a realistic vortex vacuum, where the correlations contained in vortices 
and Dirac zeromodes decay within a finite length, which should ideally be much smaller than the 
extent of the lattice. Nevertheless, a better understanding of the sensitivity of Dirac operator 
eigenmodes to Polyakov (Wilson) lines in realistic situations should be pursued, and might shed 
more light on the mechanism of chiral symmetry breaking. We plan to study this in future work. 
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